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Abstract 

Macro-orbital representation of a particle in a many body system is used to 
develop the microscopic theory of a system of interacting bosons. Each particle in 
the system represents a pair of particles moving with equal and opposite momenta 
(q, — q) at their center of mass (CM) which moves with momentum K in the lab- 
oratory frame. Below certain temperature (say, Tx), these particles assume a state 
of (q, — q) bound pairs (named as SMW pairs). The A— transition is found to be a 
consequence of inter-particle quantum correlations clubbed with zero-point repul- 
sion and inter-particle attraction; it is an onset of the order- disorder of particles in 
their (/)-space followed simultaneously by their Bose Einstein condensation as SMW 
pairs in a state oi q = qo = ir/d and K = Particles at T < Tx acquire collective 
binding which locks them at < A: >= in momentum space, < r >= A/2 in real 
space and at A(/i = 2n7r (with n = 1, 2, 3, ...) in (^— space. Consequently, the entire 
system assumes mechanical strain in inter-particle bonds and behaves like a single 
macroscopic molecule. The collective binding is identified as an energy gap be- 
tween the superfluid and normal fluid phases. The fractional density of condensed 
particles {nK=o{T)) varies smoothly from nK=o{T\) = to nft:=o(0) = 1-0. The 
A— transition represents the occurrence of twin phenomena of broken gauge symme- 
try and phase coherence. In variance with the conventional belief, the system does 
not have single particle p = condensate. In addition to the well known modes 
of collective motions (such as phonons), the superfluid also exhibits a new kind of 
quantum quasi-particle, omon, (a phononlike wave of the oscillations of momen- 
tum coordinates of the particles. Omons are sustained because the strain energy 
of inter-particle bonds is a function of relative momentum of particles. The the- 
ory explains the properties of He — II, including the origin of quantized vortices, 
critical velocities, logarithmic singularity of specific heat and related properties at 
quantitative level. It conforms to the excluded volume condition, microscopic and 
macroscopic uncertainty, and renders microscopic foundation to two fluid theory 
of Landau, \I'-theory of Ginzburg, the idea of macroscopic wave function of Lon- 
don, etc. The framework of this theory also helps in unifying the physics of widely 
different systems of interacting bosons and fermions. 
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1.0 INTRODUCTION 



Liquid '^He (LHE-4), - a typical system of interacting bosons (SIB), has been investi- 
gated extensively for its unique properties at low temperature (LT) and results are reviewed 
in several articles {e.g. [1-10] and other references cited therein). It transforms from its 
normal (N) phase {He — I) to superfiuid (S) phase {He — It) at = 2.17 K with the 
latter exhibiting several unique properties including superfluidity. One flnds that these 
properties can not be explained as the properties of a classical liquid and the same is 
true for the fact that helium remains liquid down to the lowest attainable temperature 
(T). Consequently, these properties arc, rightly, related to the manifestation of wave na- 
ture of particles at macroscopic level and LHE-4 provides unique opportunity to study 
this aspect. Naturally, different studies of a SIB are considered to be of fundamental 
importance. 

Soon after the discovery of superfluidity of He — II [11], London [12] proposed that the 
phenomenon arises fromp = condensate., -a macroscopically large fraction (np=o(T)) of 
^He atoms occupying a single particle state of momentum p — To this effect, London 
was guided by the well known conclusion that a system of non-interacting bosons (SNIB) 
should exhibit a transition identified as Bose Einstein condensation (BEC) at certain 
temperature {Tbec) [5, 13] and the LT phase should have non-zero np=o(^); -increasing 
smoothly from np=o{TBEc) = to np=o{T = 0) = 1.0 through np=o{T < Tx) < 1.0. The 
idea was initially criticized by Landau [14] because LHE-4 is not a SNIB. But the Bogoli- 
ubov's theory of weakly interacting bosons [15], concluding that the increasing strength 
of interactions simply depletes np=o{T), provided some ground. Consequently, different 
mathematical models such as Bogoliubov prescription [15], pseudo-potential technique 
[16], Jastrow formahsm [17, 18], Feenberg's perturbation [2, 18], etc. have been used 
to find the magnitude of np=o{T) in He — II and at the same time develop a general 
framework of microscopic theories hereafter known as conventional theories (CTs). While 
several articles published in [10] discuss the role of BEC to understand the LT behavior 
of widely different many body quantum systems, papers by Moroni et. al. [19] and Sokol 
[20] make useful conclusions in relation to our conventional understanding of LHE-4. It 
appears that CTs follow two different approaches, Al and A2. 

Al : This approach developed by Bogoliubov [15], Beliaev [21, 22], Hugenholtz and Pines 
[23], Lee and Yang [24], Brueckner and Sawada [25], Wu [26], Sawada [27], Gavoret and 
Nozieres [28], and Hohenberg and Martin [29], DeDominicis and Martin [30] has been 
reviewed by Woods and Cowley [4] and Toyoda [31]. While it is elegantly introduced 
by Fetter and Wadecka [32], its important aspects are summed up, recently, by Nozieres 
[33], and Huang [34]. A theoretical formulation following this approach starts with the 
hamiltonian written in terms of second quantized Schrodinger fields and proceeds by using 
important inferences such as: (i) the hard core (HC) potential can be used perturbatively 
by using a method such as pseudopotential method [16], (ii) coupling constant can be 



2 



related directly to the two body scattering length (cr) [22], (iii) a dimensionless expan- 
sion parameter (ncr^) can be used to make perturbative calculation [31], and (iv) the 
"depletion" of np=o{T) can be treated in the conventional perturbation by introducing 
c- number operator for p = bosons [35]. Using the quantum statistical mechanical ex- 
pectation value of the second quantized field operator as an order parameter (OP) of the 
transition [23, 29], it then calculates the relevant part of free energy as a function of this 
OP and explains superfluidity and related properties. To this effect, one uses sophisticated 
mathematical tools as discussed in [31, 36]. 

A2 : This approach aims at finding the radial distribution function [g{r)] and liquid 
structure factor [^'(Q)] which can be used to calculate the ground state (G-state) prop- 
erties, excitation spectrum, different thermodynamic properties and equation of state. 
Superfluidity and related properties are explained in terms of one body density matrix 
p(r) whose asymptotic value at large r gives the condensate fraction np=Q(T) [37]. While 
different aspects of this approach and related subject are discussed elegantly by Fccnbcrg 
[2], Croxton [38], Ceperley and Kalos [39] and Yang [37] and its application to LHE-4 
has been reviewed by Smith et. al. [40(a)], Campbell and Pinski [40(b)], Schmidt and 
Pandharipande [40(c)], Reatto [40(d)], and Ristig and Lam [40(e)], several methods of 
computer calculation of g{r), S{Q) and related properties are discussed and reviewed in 
[41-43]. The approach has been used by Moroni et. al. [19] and Kallio and Piilo [44] to 
determine the properties of LHE-4 and electron gas, respectively. While a comprehensive 
list of important papers related to this approach are also available in [19, 44], develop- 
ments in our conventional understanding of quantum fluids can be seen in most recent 
books and reviews [45-48]. 

Over the last seven decades a large number of papers have been published [49] to 
estabhsh the existence of np=o{T) in He — II and to flnd its value. However, several 
prominent scientists working in the field have expressed their doubts about the success of 
these approaches in concluding the desired theory. For example, reviewing the progress of 
theoretical work at the fifteenth Scottish university summer school (1974), Rickayzen [50] 
states, "There is no microscopic theory of superfluid ^ife. There are many mathematical 
models which appear to provide insight into the behavior of superfluids but there is no 
theory which provides quantitative prediction that agree with observation. Even some 
of the most widely held assumptions of the theory such as the idea of condensation in a 
zero momentum state can not be said to be proved beyond reasonable doubt." In their 
review on the subject. Woods and Cowley [4] also observe, "Despite all the experimental 
information and the numerous theoretical discussions there is still no convincing theory of 
the excitations which begins with the known interaction between helium atoms." Interest- 
ingly, though Klcban [51] proved that theories assuming the existence of p = condensate 
contradict the excluded volume condition (EVC),- a direct consequence of HC nature of 
particles but his inference did not receive due importance. Similarly, Sokol [20] makes 
interesting observations which doubt the accuracy of the inference about the existence of 
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p = condensate in He — II. He observes that: (i) neutron inelastic scattering experi- 
ments have so far not given any indication of the direct evidence {i.e. a 5-function peak in 
the momentum distribution (np) at p = 0) of the existence of np=o{T) in He — II and it is 
unhkely that this goal will ever be reached, and (ii) in the absence of detailed microscopic 
theory the estimate of np=o(0) ~ 0.1 from different experiments depends on current the- 
ories, models, and empirical relations used to interpret experimental results and if their 
underlying assumptions are incorrect, such estimates of np=o{T) will have no meaning. In 
addition though significant progress has been made in refining our mathematical methods 
and computational techniques applied to many body systems (including liquid '^He) as 
one finds from most recent books and reviews [45-48] but the presumed existence of p = 
condensate remains central to our present understanding of the superfiuidity of He — II 
and similar SIB(s); this is also evident from the most recent books [52] summarizing the 
important aspects of LHE-4. 

As such one finds that the best possible model that could be used to understand su- 
perfiuidity is the \1'— theory [53] (a union of the two fiuid model of Landau [14], quantized 
circulation [54], vanishing of superfluid density at the boundaries of the system [5], etc.) 
which has emerged as exceedingly successful phenomenology. However, the microscopic 
theory of a phenomenon has greater significance because only its agreement with experi- 
ments testifies the validity of the basic theory of nature, the above mentioned observations 
motivated us to look for an alternative approach. To this effect our initial reports [55-58] 
projecting the basic aspects of our approach greatly helped us in crystallizing our scientific 
ideas, refining our physical arguments and mathematical formulations and laying down 
the basic foundations for the final form of our theory that we report in this paper. 

Our theory is based on the macro-orbital representation of a particle which we con- 
cluded in our recent paper [59] related to the wave mechanics of two HC particles in 1-D 
box and used in our studies related to the unification of the physics of fermionic and 
bosonic systems [60], ground stated of HC core particles in 1-D box [61] and basic 
foundations of the microscopic theory of superconductivity [62]. Our theory makes no 
assumption about the basic factors {viz. the existence of p = 0-condensate or bound 
pair formation) which could, presumably, be responsible for the superfluidity of a SIB. 
In stead we only use the solutions of our A^— body Schrodinger equation to reach our 
conclusions. In fact our theory, for the first time, use the results of an obvious wave 
mechanical superposition of two particles to conclude the orderly placement of particles 
in phase space of the system and this simple effort helped us in finding an almost exact 
theory of a SIB which explains all unique properties of He — II. In variance with CTs, 
our theory concludes that the origin of superffuidity and related properties of a SIB lies 
with the BEG condensation of (q, — q) bound pairs in a state of their center of mass (CM) 
momentum K = 0. 

The paper has been arranged as follows. Analyzing the A^-particle microscopic quan- 
tum hamiltonian. Section (2.0) identifies a pair of particles as the basic unit of the system. 
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uses macro-orbital representation (Appendix -A) of a particle in writing iV— particle state 
function and concludes G-state configurations of a SIB and related aspects. While Section 
(3.0) analyzes the thermodynamic behavior of the system, Section (4.0) reports two differ- 
ent approaches used to study thermal excitations. While the origin of (q, -q) bound pair 
formation, energy gap, etc., are studied in Section (5.0), the consequences of the gap are 
examined in Section (6.0) and the consistency of our theory with: (i) phenomenological 
theories and (ii) experimental results on LHE-4 is analyzed in Section 7.0. Finally, we 
present some important concluding remarks in Section 8.0 

2.0 Basic Foundations of Our Theory 

2.1 Hamiltonian 

A system of N interacting bosons such as liquid '^He can be described, to a good 
approximation, by 

N 

H{N)^J2hi + J2V{nj) with hi^-{n^/2m)s7l (1) 

i i<j 

where m is the mass of a particle and V{rij) is a two body central force potential which 
can be expressed as the sum of: (i) a short range strong repulsion V^{rij) and (ii) a weak 
attraction V'^{rij) of slightly longer range. To a good approximation, while V^{rij) can 
be equated to hard core (HC) interaction Vnci'^^ij) [defined by VHci^ij < a) = oo and 
Vncifij > cr) = with a being the HC diameter of a particle], V'^{rij) can be replaced 
by a constant negative potential, say —Vg. Although, this seems to imply that the main 
role of V'^ij'ij) is to keep particles within volume (V) of the system but as concluded in 
Section (5.0) it also plays an important role in binding two HC bosons having equal and 
opposite (q, -q) momenta. However, to begin with, our system can be identified as an 
ensemble of HC bosons confined to a volume V and its hamiltonion can be written as, 

N N 

Ho{N) = J2h. + j:Mnj)- (2) 

i i>j 

where we use 

VHc{rij) = AS{r,j), (3) 

with A being the strength of Dirac delta repulsion. For impenetrable HC particles, we 
have A ^ oo when 0. While the equivalence expressed by Eqn.3 is mathematically 

shown by Huang [63], its physical basis can be understood by examining the possible 
configuration of two HC particle {say, PI and P2) right at the instant of their collision. 
When PI and P2 during a coUision have their individual CM located, respectively, at 
rcM(l) = cr/2 and rcM(2) = —a/2 (with tcm being the distance of the CM of a particle 
from the CM of the pair of PI and P2), they register their physical touch at r = and 
their encounter with Vffc(^ij) is a result of this touch beyond which the two can not be 
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pushed in. The process of coUision only identifies this touch; it does not register how 
far are the CM points of individual particles at this instant. In other words the rise and 
fall of the potential energy of PI and P2 during their coUision at r = is independent 
of their a and this justifies Vncif^ij) = ^^if^ij)- It may, however, be mentioned that this 
equivalence will not be valid in accounting for certain physical aspects of the system {e.g., 
the volume occupied by a given number of particles) where the real size of the particle 
assumes importance. 

2.2 Basic unit of the system 

In what follows from Section 2.1, the motion of each particle in the system can be 
described by a plane wave 

lip(b) = Aexp(ip.b) (4) 

unless it collides with other particle(s). Here A is normalization constant, while p and b, 

respectively, are momentum {wave number) and position vectors of the particle. However, 
the plane wave description is modified when the particle collides with its neighbors or 
boundary walls of the system. A collision of a particle could either be a two body collision 
or a many body collision {e.g., two mutually colliding particles also collide simultaneously 
with other particle(s)). In the former case two coUiding particles (PI and P2) simply 
exchange their momenta pi and p2 or positions bi and b2 without any difference in the 
sum of their pre- and post-collision energies. However, in the latter case PI and P2 
could be seen to jump from their state of pi and p2 to that of different momenta p'j^ 
and P2 (possibly of different energy) but it is clear that to a good approximation they 
remain in one of the possible states of two HC particles moving in the absence of other 
particle(s). Evidently, the complex dynamics of the system can be described, to a good 
approximation, in terms of the simple dynamics of a pair of HC particles as its basic unit 
and this is consistent with the fact that inter-particle forces in our system are basically 
two body forces. We analyze the wave mechanics of a pair of HC particles in Appendix-A 
and use its inferences in the following analysis. 

2.3 State function 

As concluded in Appendix-A, each particle in our system is represented more accu- 
rately by a macro-orbital {of. Eqn. A-17) rather than by 'Up(b). Accordingly, each particle 
is a member of the pair of particles moving with equal and opposite momenta (q, - q) 
at their CM which moves with moment K in the laboratory frame. It has two motions, 
q and K. While q defines its quantum size (A/ 2 = ir/q) representing the size of the real 
space occupied exclusively by it, K defines its motion as a free particle of mass 4m. By 
using N macro-orbitals for N particles, we obtain 

KiN) = nfUn) ^[(±l)^Hf expz(PK,.R,)] (5) 
p 
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as the state function which represents one of the N\ microstates of the system. Here we 
have C<ji(^i) = sin (qj.rj) and J2p' referring to the sum of A^! product terms obtainable by 
permuting N particles on different states with and (—1)''^, respectively, used for 

selecting a symmetric and anti-symmetric wave function for an exchange of two particles. 
In principle, the permutation of N particles on different states also gives N\ different 
^^(A*") and this renders 

MN)^^T.mN) (6) 

as the complete wave function of a possible quantum state of the system. 
2.4 State energy 

Since each macro-orbital representing a particle (say i—th) in $n(-^) state has a ref- 
erence to the CM coordinate system related to the pair of which the i—th particle is a 
member, we need to recast Ho{N) (Eqn. 2) in terms of these coordinates in order to 
examine if ^n{N) is an eigenf unction of Ho{N). To this effect, we define 

h(i) = lih, + h^A = ~vl.~vl (7) 

with /lAT+i = hi and write Eqn. 2 as 

N N 

^o(iv) = EMO + E^^M- (8) 

i i>j 

This easily renders [58-62] 

N 

<^rr{N)\J2A6{u^)\^^{N) >=0 (9) 

i>j 

and 

N 

En =< ^niN)\HoiN)\^niN) >= J2 

i 

which concludes that ^n{N) is an eigenstate of Ho{N) with En as its energy eigenvalue. 
While, at the first glance, Eqn. 10 indicates that En is purely a kinetic energy but the 
fact, that q values are constrained to satisfy q> qo = T^/d {cf. Appendix-A, Eqn. A-13), 
indicates that HC interaction has an obvious control on E^, of course this impact comes to 
surface only at low temperatures at which quantum effects start dominating the behavior 
of the system. 



2m 
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2.5 G-state configuration 

In the light of our inference that q > 71 /d (cf. Appendix- A, Eqn. A- 13) and K can 
have any value between and oo, the G-state energy of a SIB can be obtained by using 
all Qi = ir/di and all Ki = in Eqn. 10 and this renders 

N N N 

Eo^Y.^;;;i2=Y.l 273 ^^^h Y.'^i = V (constant) (11) 

1/3 

where we use di = with Vj being the volume of the real space,- exclusively occupied by 
the i— th particle. In writing X^f^Vj = V, we use the fact that each particle of the lowest 
possible q has largest possible quantum size A/2 and, therefore, occupy largest possible 
Vj. While different particles can, in principle, occupy different volumes but the simple 
algebra reveals that Eg is minimum value for Vi = V2 = ..vn — y/N and we have 

Eo = Nh'^/Smd'^ = Neo. (12) 

Since this means that all particles, identically, have q — Qo and K — 0, use of these values 
in Eqn. 6 renders 

$„(Ar) = nf C,o(r.) = nf sin(g„r,) (13) 

as the G-state wavefunction. It is clear from Eqns. 12 and 13 that each particle in the 
system represents a particle trapped in a box (cavity formed by its neighboring particle) 
of size d and it rests at the central point of this cavity. Since each sin(q'orj) in Eqn. 13 
represents a kind of stationary matter wave (SMW) which joins with other SMWs of 
neighboring particles at the boundaries between the two cavities they occupy, the G-state 
wave function seems to be a macroscopically large size 3-D network of SMWs which mod- 
ulates the relative positions of two particles in phase space and extends from one end of 
the container to another. In Section (5.0), we find that this network gets energetically 
stabilized due to some kind of collective binding and assumes different aspects of macro- 
scopic wave function of S-state as envisaged by London [12]. Using this information with 
our results [Eqns. A-13 and A-14 of Appendix-A)] related to < r > and < > (the 
expectation of relative position of two particles in r— and 0— spaces), we find that the 
G-state configuration of particles can be described by 

< k >= < r >= d and < cj) 2n7r with n = 1, 2, 3... (14) 

Evidently, particles in the G-state cease to have relative motion, coUisional motion and 
inter-particle scattering. Since all of them have identically equal < r >= d, (nearest 
neighbor distance) with constant distance of 2tt in the 0— space, the G-state represents 
their close packed arrangement in real space and they are constrained to move only in 
order of their locations. In the light of these inferences and our conclusion that particles 
in this state have a kind of collective binding at T < T;^ {cf. Section 5.0), it is evident 
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that particles in the LT phase have mutual binding in all the three spaces (r— , q— and 
0-). 



2.6 Evolution of the system on cooling 

In view of the condition, A/2 < d {cf. Appendix-A, Eqn. A-13), {d - A/2) for a SIB 
of fixed number density decreases with decreasing T and at certain T — Tc it vanishes at 
large. Consequently, the system at Tc has q — Qo for all particles and its state (Eqn. 6) is 
now expressed by 

= MN) E[(±l)''nf expz(PK,.R,)] (15) 
p 

For the reasons which become clear in Sections 5.0 and 6.0, the superscript S in ^^{N) 
refers to S-phase. Eqn. 15 imphes that all the A^! micro-states of our system appearing 
in $n(A^) (Eqn. 6) merge into one and the entire system at Tc attains a kind of oneness 
as envisaged by Taubes [64]. This also reveals that: (i) the system at has two separate 
components, say Fl and F2. While Fl (described by J2p' ■■ part of $^(A^)) represents a gas 
of non-interacting quasi- particle excitations originating from the plane wave AT— motions 
of particles, F2 (described by $o(A^)) represents the system in its G-state (or T = 
state) which has no motion except zero-point g— motions. The density of Fl, obviously, 
decreases with cooling the system and reaches zero value at T = 0. In Section-7.0, this 
separation is used to explain the two fluid behavior of Phase-II. 

2.7 Quantum Correlations 

Quantum correlations, basically, originating from the wave nature of particles play 
an important role in relation to the behavior of a quantum system. These correla- 
tions can be expressed in terms what is known as quantum correlation potential (QCP) 
which can be obtained [65, 66] by comparing the partition function (under the quan- 
tum limits of the system), Zq = X]nGxp(— £'„/fcBT)|$„(5')|^ and its classical equivalent, 
Zc = J2n^w{~En/kBT) exp^—Un/ksT). Here $„(S') is given by Eqn. 15. The applica- 
tion of this method to our system is justified because our theory describes the system by 
symmetrized plane waves rendering < Vhc(^) which indicates that the HC potential 
is completely screened out. Simplifying C/„, one easily finds that pairwise QCP has two 
components. The U?j pertaining to k motion controls the (p — kr position of a particle 
and we have 

C/^. = -A;BToln[2sin2(0/2)] (16) 

where T has been replaced by Tq because T equivalent of k motion energy at all T <T\ 
is To. 

U^j has minimum value {—kBToln2) at = (2n -|- l)7r and maximum value (= oo) 
at = 2mr occurring periodically at — 2mr (with n — 1,2,3, ...). Since Ufj always 



9 



increases for any small change 50 in at its minimum value, as 

11 1 

-C{S(pf = -kBTo{5(j)f with force constant C = -A;bT<, (17) 

and particles experience a force = — C50 which tries to maintain 50 = and the order 
of particles in 0-space is sustained. Since is not a real interaction like V{r) which can 
manipulate the 0-space order is, therefore, achieved when cooling drives all q towards 
qo- 

The second component pertaining to i^T— motions can be expressed by [56 and 65] 

Uij = -kBTXn [1 + exp {-2t:\E! - i?"|7A^)] (18) 



with = h/ ^2TT{Am)kBT being the thermal de Broglie wavelength pertaining to fC— motions 
for which each particle appears to have 4m mass. Note that Uij is identical to the quan- 
tum correlation potential for non-interacting particles [65]. It may be seen as the origin 
of the force that facilitates BEC of particles at the K — {) state by driving them towards 
K — Q point in space where Uij has its minimum value — /cBTln2. 

3.0 Thermodynamic Behavior 

3.1 Equation of state 

In what follows from Eqn. 10, we can express the energy of a particle in our system 

as 2 2 

e = siK) + eik) = + ^ (19) 

^ ^ ^ ^ 8m 8m ^ ' 

which can have any value between and oo. Interestingly, this possibility exists even if 

h'^k'^/Sm is replaced by the lowest energy So of its motion since K can have any value 

between and oo and we can use 

in the starting expressions of the standard theory of BEC [67, Ch. 7] to obtain 

PV 

^ = -E,(^) In [1 - zexp i-mK) + £o])] (21) 

and ^ 

'^^^^(").-exp(/3[s(K)+e„])-l ^^'^ 
with (3 = and fugacity 

z — exp {P/i) {/I — chemical potential). (23) 
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Once again, by following the steps of the standard theory of BEG [67] and redefining the 
fugacity by 



z' — z exp (— /?£o) = exp [/?(// — £o)] — exp [/?//'] with ji' — ji — So (24) 
we easily have 

P 2'K(SmkBTf'^ 



and 



N-No 27r(8m/cBT)3/2 roc ^1/2^^ i 



/oo 1 
^ In (1 - z'e--)dx = -^9,/2{z') (25) 



/oo X ' riT 



where x = (3s{K), A = /i/(27r(4m)/csT)^/^ and gn{z') has its usual expression. This 
reduces our problem to that of non- interacting bosons but with a difference. Firstly, 
we have m replaced by 4m and z by z' . Secondly, the theory of non-interacting bosons 
concludes z = 1 (or /i = 0) for T <T\ and z <1 (or < 0) for T > Tx, while our theory 
of interacting bosons fixes z' = 1 [or fi' = rendering = Eq (Eqn. 24)] for T <T\ and 
z' < 1 (or /x' < demanding 11 < So) for T > Tx. In other words we have z' and /i' in 
place oi z or II used in the theory of non- interacting bosons [67]. As such we can use 
Eqns. 21 and 22 and Eqns. 25 and 26 to evaluate different thermodynamic properties. 
For example, using Eqns. 25 and 26, we find the internal energy U — ~^{-^f)\z,v of the 
system. We have 

U = ^kBT^95/2{z') + Nso ^U' + Nso (27) 

where U' — ~-§p{-^f)\z',V: obviously, represents internal energy of motions, while 
NEo comes from /c— motions. Similarly, the Helmholtz free energy of the system can be 
expressed as 

F = Nfi- PV = NEo + {Nfx' - PV) = NEo + F' (28) 

with F' referring to that of non-interacting bosons. Following the standard methodology, 
we may now analyze F for the physical conditions for which it becomes critical and leads 
to superfiuidity. 

3.2 Onset oi K — Q condensate and Tx 

In what follows from Eqn. 28, F is the sum of: (i) F' — F{K), representing the 
contribution of if— morions which define a system of non-interacting quantum quasi- 
particles of {hosonic nature) and 4m mass and (ii) NEo = F{q) representing the zero-point 
energy of g— motions. Following the standard theory of noninteracting bosons [67], F[K) 
is expected to become critical at 

^''=i^^^^ = 8;^(2iTv)' ^^^^ 



11 



with the onset of BEC of particles in their K = Q state. Tbec in Eqn. 29 represents the 
usual point of BEC in a system of non-interacting bosons and | factor signifies that each 
boson for its motion behaves like a particle of mass 4m and Tbec varies as — . However, 
this onset would occur only when the particles at large have fallen into the G-state of 
the other component of their motion {viz. g— motions) because in an otherwise situation 
/T— motions can gain energy from g— motion energy. In what follows from Section 2.6, a 
SIB assumes a state oi q = Qo for all particles at 

h'^ 1 

which represents the T equivalent of the G-state energy Eq- It is evident that as soon as 
all particles have q = Qo state, further cooling of the system would immediately trigger 
the process of their fall in = state. Evidently, the effective temperature, at which 
particles have the state oi q = qo and the onset oi K — condensate, should be 



87Tm*kB 



1 / iV \i 
rf2 + U.61V 



(31) 



Here wc replace m by m* to account for the impact of physical parameters such as pressure 
(P) on Tx. In this context it may be noted that, in spite of the fact that particles in our 
system to a good approximation represent HC particles moving freely on a surface of 
constant —Vg, one can not ignore the role of inter-particle interaction during their relative 
motion. We note that d, V and m* are three quantities which may change with increasing 
P. While Tx is expected to increase with increasing P for the usual decrease in the values 
of d and V, however, Tx may show a reverse change if m* increases with P. In this context 
we note that m* for ^ife atoms in LHE-4 should increase with P for an obvious increase 
in the strength of inter-particle attraction with increasing P. Evidently, Eqn. 31 can 
explain the P dependence of Tx of LHE-4. 

3.3 Nature of Transition 

When our system is cooled through Tx, its particles move from their state of g > vr/rf 
(i.e., (j) > 27i) to that of g = rr/d {i.e., (j) = 2tt) {cf. Section-2.5). Evidently the system 
transforms from a state of random distribution of its particles in 0— space to that of orderly 
distribution with = 2n7r (n = 1,2,3, ..). This means that A— transition represents an 
onset of order-disorder of particles in space accompanied by the BEC of particles in the 
state of K = and q = n/d. Since BEC of particles in K = state is not different from 
the BEC of non-interacting bosons (a well known second order transition [67]) and their 
order-disorder basically represents their self-organization in 0— space which results from 
a reshuffle of their momenta, it is evident that A— transition is a second order transition. 
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3.4 Free energy and order parameter of Tx 



In the light of Section 3.1, it is evident that free energy F of the system has two 
components and it can be expressed as 



F = F{q) + F{K) ^ NEo + F{K) 



with 



F{K) = kBT 



r 

J c 



X 



z e 



(32) 



(33) 



J 

Eqn. 32 would also hold for T > Tx, provided such a T is low enough to allow only insignif- 
icantly low number of particles in the excited states of momentum q > 2qo. Guided by the 
experimental fact that a system like LHE-4 exhibit superfluidity and related properties 
below Tx, F can be expressed as 



^ 4r O 



(34) 



with Q being the order parameter, Fg being the constant component of F and A, B, C, 
etc. being the coefficients of its expansion which may depend on physical conditions such 

as T and P. However, it may be mentioned that it is F{q) (not F{K)) which is basically 
responsible for the unique properties of a SIB below Tx since a gas of non-interacting 
bosons is not expected to exhibit superfluidity and related properties in spite of its BEG 
into zero momentum state. Defining A — a.{T — Tx)/Tx, we find that 



n{T) oc 



Tx-T 



Tx 



(35) 



However, since To represents the T equivalent of the G-state energy of a particle £„, there 
is hardly any excitation in the system at Tq. Gonsequently, the state of the system at T^ 
is effectively a,T — state and for this reason we can renormahze the temperature scale 
by replacing T in Eqn. 35 by T* = T - T^ and recast Q{T) as 



n{T*) oc 



i Tx-T 
T.-T^ 



(36) 



which rightly indicates that VL{T*) reaches its maximum value 1.0 at T = T^. 



Now we may identify the physical nature of Q and to this effect our first choice is the 

fraction of the number of bosons {nK=o{T) = Nk=o{T)/N) condensed into the G-state 
defined hj q = qo and K = 0. Applying the standard theory of BEG of non-interacting 
bosons [67] to Fl, we find that nK=Q{T) should vary as [1 — (T/T^)^/^]. However, for our 
system nK=o{T) is replaced by nK=o{T*) which varies as 



nK=o{T* 



Nk=o{T*) 
N 




T-Tq 
Tx-Tr, 



3/2- 



(37) 
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Our next choices could be: (i) the strain in the interparticle bonds [{dr — d\)/dx] whose 
origin, relation with (q, — q) bound pair formation, volume expansion on cooling through 
Ta, and energy gap has been discussed at length in Section-5.0, and (ii) n* — [N*{Tx) — 
N* (T)] / N* (Tx) which represents the fraction of particles that move from their state of 
disorder in 0— space to that of order on cooling though Tx {cf. Eqn. 58 for an expression 
for N*{T)). We analyze n* and its relation with the energy gap in Section-5.0. One may 
find that all these quantities are inter-related with each other. 

3.5 Single petrticle density matrix and ODLRO 



Using Eqn. 15, we obtain the single particle density matrix. 



p(R* - R) = 



N 
V 



Nk=o{T*) , V 



N 



-27r. 



|R* -R|' 



A; 



/ 2 



2 „,7r(r"-r' 



with Nk=o{T*)/N = nK=o{T*) {cf., Eqn. 37). Here we use: (i) q<,.(r" - r')^ 
qo.(r" — r') and (ii) a renormalized T scale T* — T — Tq. 



(38) 
2mr + 



While the term in big (..) of Eqn. 38 represents the variation of density over a single 
antinodal region of the SMW of a particle, Nk=o{T*) stands for the number of particles 
condensed to the state oi K — and q — Qo', X't — h/ ^27r{Am)kBT represents thermal 
wavelength attributed to K motions. Wc note that under the limit |R* — R| ^ oo. the 
"one particle density matrix" (p(R* — R)) has nonzero value {Nk=o{T)/V) for T < Tx 
and zero value for T > Tx since Nk^o{T*) is 7^ for T < Tx and = for T > Tx- Evidently, 
our theory satisfies the criterion of Penrose and Onsager [35] for the occurrence of BEG 
in the G-state of the system defined hy K — and q — n/d and agrees with the idea 
of ODLRO, spontaneous symmetry breaking and phase coherence advanced, respectively, 
by Yang [37] , Goldstone [68] and Anderson [69] . 

3.6 Logcirithmic singularity of specific heat 

As shown in this section, this singularity is associated with the process which reshuffles 
particles from their random positions in 0— space at T^ (just above Tx) to their ordered 
locations at T;^(just below Tx). We assume that in this process about A^^ particles move 
from their = (2n7r ± S(f)x) positions at T^ to (f) — {2n + l)7r at T^ . This releases Ae 
energy responsible for the singularity and we have 



As = -NxkeTo 



(ln2)sin^ 



2n7r± 



In 2 



Following the theories [70] of critical phenomenon we may define 

5(t>x^5(t>x{o)\C\'\l + a2\C? + a,\Cf 



(39) 



(40) 
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with ( — '^Tx^ ■ ^ good approximation we have 

A..-.V(I^)..T„.„(^£^)^ (41) 

by using S(l)x = 50A(o)|Cr ^^'^ -^a = ^ ^t^^ ' latter expression is so chosen to ensure 
that As does not diverge at Tx and it decreases with decreasing T through 7^. Eqn. 41 
gives 

Cp(r ^ Tx) ^ -^kBTo[2u\n |C| + In (50a(o)^) - ln4 + 2u] (42) 

It may be mentioned that 54>x{o), etc. may not have same values at T > and at 
T <Tx, obviously, because the relative configuration of the system on the two sides of Tx 
differ significantly. 

4.0 Thermal Excitations. 

4.1 Feynman's approach 

Defining an excited state of the system hy = J2if{'r'i)4> and the G-state by 0, 
Feynman [71] showed that the excited state energy is minimum for /(n) = exp(ik.rj) and 
he obtained 

with S{Q) = structure factor of the system. However, E{Q) obtained from Eqn. 43 was 
found to be about two times the experimental value. Introducing back fiow effect Feynman 
and Cohen [72] later found better agreement with experiments but with considerable 
discrepancy at higher Q. In this section we use Feynman's approach to our particles in 
(q, -q) configuration. We note that under the impact of an excitation a SMW pair in 
the G-state configuration of (qo, — qo) and K = 0, expressed by (ilo) (Eqn. A-17, with 
k — 2qo and K — 0) moves to a new configuration (qo + Aq, — qo + Aq) described by 

CJr) = C(go)(r) exp (zQ.R) exp {-ie{Q)t/h (44) 

with 2Aq = Q which means that the impulse changes in the CM momentum and energy 
of both particles and its impact should be identified by Q not by Aq. Eqn. 44 further 
shows that f{R) = exp (iQ.R) is the real form of / which renders = I]i exp (zQ.Rj)0. 
Using these facts and recasting the relation, — (/i^/2m)V^/(-R) = £ J g{R — R')f{R')d^R' 
(Eqn.(11.25) of Feynman [71], p 329-30), for i-th particle we find -{h^ /8m)V'if{R) = 
£i{Q)S{Q)f{R); note that our single particle energy operator is — (^^/8m)V^. This ren- 
ders {h'^Q'^ /8m) — ei{Q)S{Q). Adding this relation to a similar relation for j-th particle 
of the pair gives 



h'^Q^ ^ 1 

AmSiQ) ~ 2 



E{Q) = SiiQ) + s^iQ) = = i;E{Q)Feyn (45) 
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which naturally explains E{Q)expt of He — II. 

4.2 Aproach of mono-atomic chain 

We note that a system like liquid '^He is expected to exhibit: (i) no transverse mode 
because the shear forces between its particles are negligibly small, and (ii) only one branch 
of longitudinal mode because the system is isotropic. Further since the particles in the 

G-state of a SIB from a kind of close packed arrangement with A0 = kr = 2mi (Section 
2.5), we can visualize waves of 0-oscillations. Using a linear chain of atoms with nearest 
neighbor forces, we find that the frequency dispersion of 0-oscillations can be expressed 

by 

= V(4C)//3|sin(Qd/2)| (46) 

where Q is the wave vector and j3 is the measure of inertia for motion. Since 0— positions 
of particles are restored by C/?,, C is given by Eqn. 17. We note that 0-oscillations can 
appear as the oscillations of r and q because 5(j) = 2qSr + 2Sqr; we have phonons when q 
remains unchanged at qo, and omons (a new kind of quantum quasi-particle representing 
a phononlike wave of the oscillations of momentum) coordinates when r remains fixed at 
d. 

Evidently, uJr{Q) of phonons can be represented, to a good approximation, by the 
dispersion of the clastic waves in a chain of identical atoms and it can be obtained from 
Eqn. 46 by replacing (3 and C by m and C*, respectively. We have 

C* = An^C/d'' = 27r^kBTjd^ = i^K'/Amd'' (47) 

However, for a better accuracy d and C* should, respectively, be considered descending 
and ascending functions of Q because increase in the energy of particles changed by an 
excitation reduces their WP size and decreases d and thereby increases C. Consequently, 
EphiQ) can be expressed more accurately by 

EphiQ) = huJriQ) = h^AC{Q)/m\ sm{Qd{Q)/2)\ (48). 

We note that use of C*{Q) and d{Q) not only explains the experimentally observed Eph{Q) 
of He — II but also accounts for its anomalous nature at low Q [56]. This aspect has 
been studied in detail in [73, 74]. However, since d{Q) can not be smaller than a, d{Q) 
and C{Q) are bound to become Q independent for Q > t:/(t and the maximum in Eph{Q) 
{i.e. the position of so called maxon) should fall at maxon Q — Qmax = tt/it and Eph{Q) 
ior Q > n/a and < 27r/d should follow 

EphiQ) = ^^riQ) = n^^C{Qraax)/m\ sm{Qa /2)\ (49) 

It is evident that phononlike dispersion is expected till the excitation wavelength A remain 
larger than d {i.e. Q < 2T:/d). However, the momentum and energy of the excitation 
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would be carried by only a single particle, if Q < ^n/d {i.e. A < cr) and this implies that 
E{Q) for Q >2Ti/d would follow 

Esp{Q) - n'Q''/2mF (50) 

which represents a kind of single particle dispersion with mp being a kind of effective 
mass that measures the effect of quantum correlation in (q, -q) pair configuration. It is 
expected to be around 4m near Q — 2Tr/d and then decrease slowly with increasing Q 
beyond this point. 

The transition of E{Q) from Eph{Q) to E.^p{Q) would, obviously, occur at 27r/d < Q < 
2ir/a, before the Eph{Q) (Eqn. 49) meets its zero value at 27i/a. This implies that E{Q) 
has to have its minimum {identified as the roton minimum for He — II) at a Q = Qmin 
near the mid-point oi Q = 2'n/d and Q = 27i/a, i.e. 

(51) 

Evidently, Eqns. 48-51 represent a Landau type spectrum and use of Eqn. 47 in Eqn. 
46 renders 

Vp — Vg — \fKh/2md (52) 

where Vp and Vg, respectively, represent the phase and group velocities of phonons at 
g fa 0. 

We further find that the equation of motion of rg (the r of s-th atom) 

d^rs = -^ujl[2rs - rs-i - r^+i] (53) 
transforms into a similar equation for ps = ^<ls by operating mdt and this concludes 

Ur{Q)^Ug{Q) (54) 

implying that the omon dispersion u}q{Q) is not different from phonon dispersion u}r{Q). 
Interestingly, as concluded in Section (5.0), an omon is an anti-phonon quantum quasi- 
particle. 

5.0 Energy Gap, Self Energy and Bound Pairs 

The system ai T < T\ is separated into two components, Fl and F2 (c/. Section 
2.6). While Fl depicts i^— motions as quasi-particle excitations whose density decreases 
on coohng and vanishes at T = 0, F2 represents the G-state of the system and portrays 
zero-point g— motions with q — q^. Evidently, the unique low T properties of our system 
that last even at T = 0, are related to F2 and to this effect we examine its evolution on 
cooling the system from T;^ to T = 0. 



71 n 
d a. 
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As concluded in Sections 2.5 and 2.6, particles in the system are locked at eo{Tx). 
However, when the system is cooled to T < Tx, = h'^/Smd'^ tends to have a value 
lower than £o(2a) which, obviously, requires an increase d from dx to dr — dx + Ad; 
here dx means d at A— point. Naturally, this brings zero-point force (/„ = —ddSo ~ 
K^/Amd^) into operation which increases d by pushing two neighboring particles to an 
increased distance dx against their ' inherent attractive force /„. Consequently, the system 
has volume expansion on cooling around Tx and this agrees exactly with experimentally 
observed negative thermal expansion co-efficient of hquid ^He a,t T < Tx [1] . Evidently, 
the evolution of F2 on cooling is an intcr-play of /„ and fa- In the following wc Tise 
this inference and other important aspects of F2 to present three inter-related pictures to 
conclude the formation of (q, - q) bound pairs, energy gap and self energy. 

(1) : The impact of /„ on the states of the system at T < Tx, which is decided by its balance 
with fa, can be understood by using V{rij) as a perturbation. We first demonstrate this 
for the states of a pair of particles and to this effect diagonalise (2x2) energy matrix 
defined by En = E22 = £0 and E12 = E21 = Po with Po =< ^(^ij) >• We note that fo is a 
consequence of Vhc'('") clubbed with wave nature of particles which increases its operation 
from cr to A/2, particularly, for the particles of A/2 > a [59]. Naturally, it appears on the 
scene only at T < Tx where A/2 ^ rf^ (little more than d decided by V{rij). The macro- 
orbitals of particles under these physical situations should, obviously, have their overlap to 
render non-zero < V{rij) > and this is in line with our zero-order results < Vncif^) 0- 
Evidently, one finds that the pair has two states of energy ± in fact \Po\ should 
better be rei)laced by |/9o(T)| as the overlap of the macro-orbitals of two particles may 
depend on T. The states of energy {Eo — |/9o(T)|) and (so + |/?o(T)|) can, respectively, be 
identified as bonding (or paired) and antibonding (or unpaired) states. This follows the 
established approach of Molecular Orbital Theory [75] applied to a similar case in which 
two identical atomic orbitals form two molecular orbitals of bonding and anti-bonding 
nature. The pair is expected to be in bonding state provided the two particles remain 
locked in the relative configuration characterized by Eqn. 14. 

Now we apply the same approach to the state of N particles by constructing a A'^xA'" 
matrix for Hr{N) (representing the relative motion of A" particles) with [Hr{N)]mn — £0 
for m = n, and [Hr{N)]mn = [y{i^)]mn ^or m ^ n with [K(r)]m„ having non-zero value 
only if m and n refer to two neighboring particles. Note that each particle has 6-12 
(depending on the symmetry of their assumed spatial arrangement) nearest neighbors. 
The diagonahsation of this matrix renders N/2 energy levels with energy > Neo{Tx) {the 
anti-bonding states) and N/2 energy levels with energy < Neo{Tx) {the bonding states). It 
is natural that F2, separated from Fl, assumes the lowest possible bonding state. Since 
the perturbativc effect of V{r) lowers Eo of each particle, identically, all particle fall in the 
bonding state simultaneously and acquire a kind of collective binding. 

(2) . Although, the preceding analysis renders a good account for the origin of the col- 
lective binding of particles, however, a relation that helps in calculating its magnitude is 
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needed. To this effect, we note that the energy of F2 deceases from Neo{Tx) to N.SoiT) 
when the system moves from its A— point to T < T\. This concludes that energy of each 
particle falls by 



Ae = [eo{Tx) - £o(r)] 



dx — d\ 



(55) 



Amd\ d\ 



d\ 



A simple analysis of the equilibrium between fo and fa reveals that half of the Ae is 
stored in the system as strain energy in inter-particle bonds, while the remaining half is 
lost to surrounding. Since this happens to all particles, the net energy loss (= NAe/2) by 
F2 turns out to be their collective binding. Using the coherence property of the system 
(evident from A0 = 2nn), one may find that the effective binding per particle becomes 
much larger than KbT [76] , even when the real binding per particle is very small and this 
ensures the stability of our system with collective binding among its particles. Evidently, 
the entire system at T < seems to represent a macroscopically large single molecule 
as envisaged by Foot and Steane [77] for the BEC state of trapped dilute gases. It is 
obvious that the stability of this state can not be disturbed by any perturbation of energy 
< NAe/2 and we may define 



as an energy gap between the S- and N-states. 

The strain energy AVg = NAe/2, stored with the elongated bonds between nearest 
neighbors, obviously, represents the net increase in the potential energy (= AVg) of parti- 
cles and we call it as the self energy of the system. Since the strain {Ad/d\ = {dT — d\)/d\ 
arising due to increase in quantum size of each particle from d\ to dr) is an obvious 
function of the momentum q of each particle, we have AVs = AVs{qi,q2, ■■■qN) and this 
prepares the system to sustain phononlike waves of collective oscillations of momentum 
coordinates of particles (named as omons, cf. Section-4.2). Evidently, AVg can also be 
recognized as the energy of omon field. The fact that AVg increases with decreasing T 
implies that omon field intensity increases when phonon field intensity decreases and vice 
versa and this means that either the omon field intensity is the intensity of phonons 
condensed with the system or an omon is an anti-phonon quantum quasi-particle. We 
note that AVs assumes it maximum value at T = and serves as the source of collective 
motions at T = when phonon cease to exist. 

Since each particle in our system represents a (q, - q) pair, the above analysis concludes 
that all particles at T < Tx assume the states of (q, - q) bound pairs. Further in view 
of our results (c/., Eqn. 14), which indicate that each pair of particles in F2 is locked at 
< k >— 0, < r >— d and < (j) >— 2n7r, our bound pair represents two particles bound in 
all the three {k—, r— and 0) spaces. However, this does not imply that two particles in 
the system form a He2 type diatomic molecule; the entire system assumes a state where 



Eg{T) = NAe/2 = Ne, 



dr - dx 



(56) 
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particles interact with their neighbors identically by two body forces. We use the word 
bound pair because each particle is in a state of (q,-q) pair. 



(3). As inferred by Ulhenback and Gropper [66], a quantum system can be treated like 
a classical system by adding quantum correlation potential to the hamiltonian of the 
system. Using this inference and Eqn. 16, we find an alternative relation for Eg{T). 
We note that F2, for its proximity with Fl, has a small number of particles N*{T) in 
thermally excited g— motion states of g > 2qo. These particles are, obviously, devoid of 
quantum correlations with particles in the G-statc {q = go). However, in the process of 
cooling, these particles fall to the G-state and establish quantum correlations with other 
particles of this state. Consequently, the energy of each particle falls below zero energy 
level by —/csTo In 2. Assuming that number of these particles decreases from N*{Tx) to 
N*{T), we find 

Ae(r) = -ksTo In 2[N*{Tx) - N*{T)] (57) 
as the net fall in energy due to quantum correlations. Here we have 

- Eo 



2m 1 
7? 



exp 



knT 



1 -1 

^/ede (58) 



where Ec — h^Ql/2m (with Qc ~ 27r/(7). Once again it is evident that half of the 
Ae(T) energy released out is stored back in the system as strain energy in expanding 
nearest neighbor intcr-particlc bonds which indicates that the net fall in energy per particle 
is ^Ae(T). Since the negative value of |Ae(T) means a kind of collective binding of 
particles and we have Eg{T) = ^Ae{T) is an alternative relation for Eg{T). Its accuracy 
is corroborated by the fact that ^Ae{T) (Eqn. 57) and Eg{T) (Eqn. 56) have closely 
equal value for He — II at all T <T\ [56]. Note that Eqn. 58 is an approximate relation 
for N*{T), since it uses a free particle dispersion, e — h^Q'^/2m which is valid only to a 
good approximation. 

6.0 Energy Gap and its Consequences 

In what follows from the Section-5.0, the energy of F2 at T < Tx can be expressed 
as F{q) — Neo{Tx) — Eg{T) where NSoiTx) is constant. Evidently, the origin of different 
properties (including superfluidity and related aspects) of our system lies with Eg{T) and 
this fact is used for the following analysis. 

1. Superfluidity and Related Properties : 

If two heads X and Y in the system have small T and P (pressure) differences, the 
equation of state is Eg{X) = Eg{Y)+SAT-VAP. Using Eg{X) = Eg{Y) for equilibrium, 
we get 

SAT = YAP (59) 
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This reveals that : (i) the system should exhibit thermo-mechanical and mechano- caloric 
effects, and (ii) the measurement of rj by capillary flow method performed under the 
condition AT = and of thermal conductivity (©) determined under AP = should 
reveal r] = and 6 ~ oo, respectively. As such the S-phase is expected to be a superfluid 
of infinitely high ©. 

Interestingly, several important aspects of our system can also be followed qualitatively 
from the configuration of F2. For example, we note that: (i) a close packed arrangement 
of particles in a fluid like system can have no vacant site, particularly, because two neigh- 
boring particles experience zero point repulsion which, naturally, means that the system 
should have very large 0, (ii) the system can not have thermal convection currents for 
its large © and close packing of particles and this explains why He — II does not boil 
like He — I, (iii) since particles in F2 can move only in order of their locations and cease 
to have relative motion, the system is bound to exhibit vanishingly small r], particularly, 
for their flow in narrow capillary, etc.) In the rotating fluid, however, particles moving 
on the neighboring concentric circular paths of quantized vortices have relative velocity 
as a source of natural viscous behavior. This explains both viscosity and rotation para- 
doxes [5]. As such the loss of viscosity in linear motion is not due to any loss of viscous 
forces among the particles, rather it is the property of the S-phase configuration in which 
particles cease to have relative or collision motion. 

2. Critical Velocities and Stability of S-phase : 

Using the same argument, which renders Eqn. 44, we find that the state function 
^o{N) of the S-phase changes to ^1{N) when the system is made to flow with velocity 
Vf — hAq/m. We have 

$:(7V) = $„(7V) exp (zK.Ef Ri) exp [-i[N{eo + e{K)) - Eg{T)]t/h] (60) 

with 2Aq = K. This reveals that the S-state function remains stable against the flow 
unless its energy Nmv'j/2 = Ne(K) overtakes the collective binding Eg(T). We use 
this observation to explain critical velocity Vc for which the system loses superfluidity. 
Equating Eg{T) and flow energy {N.mvj/2 with Vf = Vc), we obtain the upper bound of 
Vc- We have 

v,(T) = ^[2Eg{T)/Nm] (61) 

Avc < Vc{T), at which the superfluid may show signs of viscous behavior, can be expected 
due to creation of quantized vortices. However, this cause would not destroy superfluidity 
unless energy of all vortices produced in the system exceeds Eg{T). 

3. Coherence Length : 

Since the main factor responsible for the coherence of F2 is its configuration which 
locks the particles at A0 = 2m: {cf. Eqn. 14) with collective binding Eg{T), the coherence 
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length (not to be confused with heahng length [5]), can be obtained from 

= l/mv,(T) = h^[N/2mE,{T)] (62). 

4- Superfluid Density : 

Correlating the superfluid density, ps, as the order parameter of the transition, with 
Eg{T) we find a new relation 

to determine Ps{T) and normal density, Pn(T) = p(T) — PsiT). Evidently, VciT), C(T), 
and ps{T) can be obtained if we know Eg{T) (Eqn. 56). Further since S-state function 
vanishes at the boundaries of the system, it is natural that Eg{T) and Ps{T) also vanish 
there. 

5. Superfluid Velocity : 

Concentrating only on the time independent part, Eqn. 60 can be arranged as 

$:(7V) = ^o{N) exp {iS{R) (64) 

with 

S{R) = K.(Ef Ri) (65) 
being the phase of the S-state. This renders 

. Av„,S(«) ^ ^ (66) 

as a relation for the superfluid velocity; here we use the fact that \/ji^S{R) renders the 
momentum of the pair (not of a single particle). One may find that Eqn. 66 does not 
differ from the superfiuid wave function presumed in the ^—theory of superfiuidity [78] 
(c/. Section (2.3) of Reference [8]), of course for the well defined phenomenological reasons, 
^—theory assumes S{R) to be a complex quantity. 

6. Quantized Vortices 

Using the symmetry property of a state of bosonic system, Feynman [54, 71] showed 
that the circulation, k, of the velocity field should be quantized, i.e., k = ^ with n — 
1,2,3,... However, Wilks [1] has rightly pointed out that this account does not explain 
the fact that He — I to which Feynman's argument applies equally well, does not exhibit 
quantized vortices. Using Eqn. 66, we find that 

K = Evs(^)•Ar. = - E Aq.-Ar. = — (67) 
, m, m 
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by using the condition that Aqj.Arj = 2n7r which presumes that particles moving on 
a closed path maintain phase correlation. To this effect our theory reveals that particles 
have their 0-positions locked at = 2mr only in S-phase which indicates that only 
this phase can exhibit quantized vortices. However, since particles of N-phase have ran- 
dom distribution (A0 > 2n7r) in 0-space, this phase can not sustain 0— correlation and 
quantized vortices. 

7. S-state and its Similarity with Lasers 

We note that the system below Tx defines a 3-D network of SMWs extending from its 
one end to another without any discontinuity. In lasers too these are the standing waves 
of electromagnetic field that modulate the probability of finding a photon at a chosen 
phase point. The basic difference between the two lies in the number of bosons in a single 
anti-nodal region (AR) of a SMW. In case of lasers this could be any number since photons 
are non-interacting particles but for a SIB like *^ife or ^^Rb one AR can have only one 
atom. 

7. Consistency with Phenomenological Theories and Experiments 

1. Two fluid theory 

Since our system at T < Tx has two separated components Fl and F2, it can be 
identified to be a homogeneous mixture of two fiuids. While Fl, described by plane 
waves of momentum K, represents a gas of non-interacting quasi-particle excitations, F2 
representing zero-point q'— motions describes the G-state of the system where particles 
are locked with < k >= 0, < r >= A/2 and Acf) — 2mr and consequently cease to 
have relative motion, or collisional motion essential for non-zero viscosity. Naturally, the 
thermodynamic properties, such as specific heat, non-zero entropy, etc., as well as non- 
zero viscosity are contributed totally by Fl. As such Fl has all properties of a N-fluid, 
while F2 having zero entropy, zero viscosity, etc. has the properties of S-fluid. Most 
interestingly our theory provides microscopic origin to the two fluid theory of Landau 
[14] and in spite of its separation into Fl and F2, each and every particle can be seen to 
participate in both fiuids. 

Theory 

We find that superfiuidity is basically a property of the G-state (F2 component) of 
our system represented by Eqn. 13 which can also be expressed as 

^o{N) = (68) 

with n — N/Y. However, when F2 is made to flow its state is given by $*(A^) = 
^o{N) QW^S{R) (Eqn. 64). For the phenomenological reasons (such as the superfluid 
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density is not always equal to n), we assume that the phase S{R) (Eqn. 65) is a 
complex quatity given by 

S{R) = ^R) + iCi{R) (69) 

which renders Ug — nexp (— 2^j(i?)). Evidently, $*(-/V) has the structure of ^—function 
that forms the basis of the well known ^—theory of superfluidity. This shows that our 
theory provides microscopic foundation to the highly successful ^f— theory [78]. 

3. Properties of He-II 

(i) Thermodynamic properties : An estimate of T\ for LHE-4 by using Eqn. 31 (with 
m* = m) renders a value falling around 2.1 K which agrees closely with experimental 
value 2.17 K. In view of an observation by Woods and Cowley [4], Feynman's relation 
E{Q)Feyn — h'^Q^/2mS{Q) [71] renders E{Q) values that are nearly two times of E{Q)expt 
(experimental E{Q) for He — II). Evidently, our E{Q) (Eqn. 45) equaling one half of 
the E{Q)Feyn, cusures good matching with E{Q),.xpt) for He-II and this implies that our 
theory can explain the thermodynamic properties of He — II accurately. In addition, 
using our alternate set of relations for E{Q) (Eqns. 46-51), concluded by our theory, we 
find good agreement between theory and experiment for : (i) for the anomalous nature 
of phonon dispersion at low Q, velocity of sound at low Q [56 and 80]. 

The problem of explaining the experimentally observed logarithmic singularity in 
Cp{T) of LHE-4 at Tx has been a challenging task (c/. a remark by Fcynman in his 
book [71, p.34]). Naturally, the fact that for the first time our theory succeeded in ex- 
plaining it at quantitative level specks for the accuracy of the microscopic foundation of 
our theory. Using the parameters of liquid ^He in Eqn. 42 with v — 0.55 and 5(f)\{o) — tt, 
we find 

Cp(J/mole.K) « -5.71 In |C| - 10.35 = -A. In |C| + S (70) 

which can be compared with experimental results Cp = —5.355 In |C| — 7.77 for T > Tx 
and Cp = — 5.1 In ICI -1-15.52 for T < Ta [79]. The fact that our A value agrees closely with 
experiments speaks of the accuracy of our theoretical result. With respect to our choice 
of = 0.55 and S(j)x{o) = tt, we note that: (i) S(f)x originates basically from change in 
momentum AA; ^ and ^ varies around Tx as |T — Txl^"; note that critical exponent i' 
lies in the range 0.55 to 0.7 [70]. For our choice of S4>x{o) = n, we note that tt is the largest 
possible value by which phase position of a particle changes. We note that experimental 
value of S = -7.77 (for T > Tx) and + 15.52 (for T < Tx) can be obtained if 5(/)a(o) is 
chosen to be around 0.787r and O.OStt, respectively. These arc consistent with the reasons 
of assymetry of this singularity as stated in Section-3.6 and the fact that 0— fluctuations 
at T > Tx are expected to be much higher than those at T < T^. 

(i) Hydrodynamic properties : Our theory provides microscopic foundations for: (i) two 
fluid theory of Landau (Section 7.1), (ii) the \t'— theory of superfluidity (Section 7.2), (iii) 
the quantized vortices (Section 6.6), (iv) the vanishing of Ps{T) at the boundaries of the 
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system (Section 6.4), (v) the collective binding or energy gap (Section 5.0) which not 
only ensures the stabihty of S-state but also accounts for the critical velocities, coherence 
length, etc. (Section 6.0). Evidently, it can explain the hydrodynamic properties of 
superfluid SIB like He — II. While a brief analysis of the agreement between theory and 
experiments for He — II is available in [56], a detailed quantitative study, being completed 
by our group [80], would be published soon. 

As shown Fig.(l), we note that the T dependence of: (i) the fraction of particles con- 
densed inK = state, nK=o(T*) [Curve A*, Eqn. 37] and (ii) n*(T) = [l-N*{T)/N*{Tx)] 
[Curve B, Eqn. 58] have an excellent agreement with similar nature of the order param- 
eter {fl{T*) [Curve C*, Eqn. 36] of a second order transition. In addition they have 
similar agreement with the T variation of : (i) ps{T) (Curve E2) deduced from the ex- 
perimentally observed second sound velocity of He — II (taken from [5]), and (ii) the 
normalized inter-particle bond strain dx — d\/do — d\ (curve El) obtained from experi- 
mental values of p{T) [1]. These facts, evidently, identify inter-particle bond strain as a 
basic aspect of the A— transition, T* as the effective temperature scale for the order pa- 
rameter, and establish the accuracy of our theory in accounting for the unique properties 
of He — II. It is further corroborated by the fact that El and E2 curves do not agree 
with nK=o{T) = 1 - {T /Txf/^ (Curve A) and Q{T) [Curve C, Eqn. 35]. As such an 
overall agreement between experimental and theoretical results for LHE-4 establishes the 
accuracy of our theory. 

8.0 Concluding Renicirks 

This paper presents an alternative approach to the microscopic understanding of a 
SIB. We find that each particle in the system has two motions, q— and K—, because it 
represents a pair of particles moving with equal and opposite momenta (q, -q) at their 
CM which moves with momentum K. Consequently, it should be described by a macro- 
orbital. The onset of A— transition is an order-disorder of particles in 0— space followed 
simultaneously by their BEC in the state oi q = Qo and K — 0. It is a consequence of 
quantum correlations which drive 5— values towards q — qo and values towards K — Q 
and an inter-play of zero-point force fo and inter-particle attraction /„ which leaves a kind 
of mechanical strain in the inter-particle bonds as the basic component of order parameter 
(Cf. Fig.l curve El). The system is expected to exhibit —ve thermal expansion coefficient 
around Tx as evinced by LHE-4. 

The transition separates the two motions and the system behaves like a homogeneous 
mixture of two fluids: (i) Fl, -representing a gas of quasi-particle excitations describing 
the plane wave i^— motions of N particles and (ii) F2, -representing the G-state of the 
system where particles are locked at < r d, = 2m: and zero-point motions at 
q — qo and K — The condensate fraction nK=Q{T) rises smoothly from nK=Q{T) — 
at T = Tx, reaches nK=o{T) ^ 1.0 around T ^ To and nK=o{T) = 1.0 at T = 0. The 
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HC interaction leading to excluded volume condition [51] pushes all particles to occupy 
identically equal volume and a state of identically equal q = go- 
Superfluidity and related properties are, basically, the properties of T = state of 
F2 component of a SIB which implies that A— transition is a kind of quantum transition 
which occurs at non-zero T for the proximity of F2 with Fl. Particles in F2 represent 
(q, -q) bound pairs, not only in gr— space but also in r— and 0— spaces. They assume a 
kind of close packed arrangement with collective binding which leads the entire system to 
behave like a macroscopically large single molecule. The collective binding represents an 
energy gap between the S- and N-fluid phases. The S-state is consistent with microscopic 
uncertainty as evident from q > re/ d as well as macroscopic uncertainty since $o(^) (Eqn. 
13) vanishes at the boundaries of the system. The difference between the S- and N-fluid 
states can be identifled with a difference in the ordered positions and motions of soldiers 
in an organized army platoon and a crowed. Our theory makes no assumption about the 
nature of BEG in a SIB. It reaches its all conclusions from the analysis of the solutions 
of its iV— body Schrodinger equations. 

In many respect, above summarized inferences of our theory differ from conventional 
theories which presume the existence of p = condensate in the S-state of a SIB. Using 
two different theoretical approaches (Al and A2 summarized in Section-1.0) and different 
possible mathematical tools, they conclude a maximum of ^ 13% p — condensate in 
He — II and 60% in recently discovered EEC state of dilute trapped gases [81] and for 
the last several decades this conclusion has been regarded as basic origin of the unique 
properties of every superfluid SIB. 

In summary, this paper provides an alternative theoretical framework concluding a 
different picture of S-phasc of a SIB like LHE-4. Analyzing the wealth of experimental 
observations on the S-state of different SIB(s), one may agree that the wave nature of 
particles has wonderful capacity for their self organization in phase space at = 2n7r 
with < r >— X/2 — d and our theoretical framework can reveal the truth of the low T 
behavior of widely different many body systems. 

Finally, as discussed briefly in [57, 60], the framework of our theory unifles the physics 
of widely different many body systems of interacting bosons and fermions including low 
dimensional systems, atomic nucleus, newly discovered BEC states [81], etc. This is 
corroborated by our detailed studies of N HC particles in 1-D box [61] and microscopic 
theory of superconductivity [62] which not only accounts for the highest Tc that we know 
to-day but also reveals that the phenomenon may, in principle, occur at room temperature 
and co-exist with ferro-magnetism. For the flrst time our theory concludes [60, 62] that 
superfluid transition in liquid ^He should occur around 2mK which matches almost exactly 
with experimental result. It may be emphasized that (q, -q) bound pairs of particles as 
concluded in [62] and the present paper should not be confused with Cooper pairs [82] 
since the two have several similarities and dissimilarities. 
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Appendix -A 

Wave Mechanics of Two HC Particles and Macro-orbitals 

detailed analysis on this topic avaible in [62] 

The dynamics of two HC impenetrable particles PI and P2 can be described by 

[-^ E V.^ + Vncir)^ 2) = E{2)^{1, 2) {A - 1) 

As argued in Section 2.1, we may use Vhc{'>^) = Ad{r) (Eqn. 3) and simplify the 
solution by using CM coordinates 

r = b2 — bi and k = p2 — pi = 2q, {A — 2) 

where r and k, respectively, represent the relative position and relative momentum of PI 
and P2, and 

R=(b2 + bi)/2 and K = p2 + Pi, (^-3) 

where R and K, similarly, refer to the position and momentum of their CM. Without loss 
of generality, Eqns. A-2 and A-3 also render 

Pi = -q+y and p2 = q+y. (^-4) 
We use these relations to express Eqn. A-l as 

(-^ Vl -- +AS{r)) *(r, R) ^ E(2)*(r, R) {A - 5) 

\ 4m m J 

with 

^(r, R) = ^fc(r) exp(iK.R). {A - 6) 

The HC interaction affects only ^/'^(r) (the state of relative motion) which represents a 
solution of ^ 

^-^ V' +A6{r)^ Mr) = E^Mr) (A - 7) 

with Ek — E{2) — h^K^/Am, while the CM motion [exp(iK.R)] remains unaffected. 

To find ^(r, i?) (Eqn. A-6), we treat A5{r) as a step potential. Since A5{r) = for 
r 7^ 0, PI and P2 can be represented, by independent plane waves but at r = where 
A5{r) = oo, we use the condition \E'(r, R)\r=o — 0. Following these points, a state of PI 
and P2 can, in principle, be expressed by 

^(1,2)^ = ^K(ri)«p,(r2) ±Mp,(ri)Mp,(r2]. {A - 8) 
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However, ^'(1, 2) (of +ve symmetry for the exchange of two particles) has to be excluded 
as undesirable function because it does not satisfy ^(r,R)\r=o = 0. But this leaves only 
^(1, 2)~ of —ve symmetry which does not fit with the bosonic character of our system. We 
encountered this problem in our recent study of the 1-D analogue of Eqn. A-7 in relation 
to the wave mechanics of two HC impenetrable particles in 1-D box [59]. Following this 
study we find that the state of PI and P2 can be expressed by 

C(r, Rf = Ck{r)^ exp (iK.R) {A - 9) 

with 

(^(r)- = v^sin (k.r/2) {A - 10) 

of -ve symmetry, and 

a(r)+ = v^sin(|k.r|/2) (^-11) 

of -|-ve symmetry. 

Cfc(r)^ represents a kind of stationary matter wave (SMW) which modulates the prob- 
ability |Cfc(^)^P of finding two particles at their relative phase (0) position = k.r in the 
0— space. Interestingly, the equality |Cfe(^") P = |Cfc(^")^P renders an important fact that 
the relative configuration and relative dynamics of two HC particles is independent of 
their fermionic or bosonic nature and the requirement of a bosonic or fermionic symmetry 
should be enforced on the wave functions of their K— motions or spin motions. We use 
this fact in constructing iV— particle wavefunction in Section 2.3. In agreement with Eqn. 
A-4, the SMW character of Cfc(^)^ reveals that: (i) two HC particles have equal and op- 
posite momenta (q,-q) in the frame attached to their CM which moves with momentum 
K in the laboratory frame, and (ii) their relative motion maintains a center of symmetry 
at their CM. This implies 

rcM(l) = -rcM(2) = ^ and kcM(l) = -kcM(2) = q {A- 12) 

where rcM(i) and kcAf(i), respectively, refer to the position and momentum of i—th 
particle with respect to the CM of two particles. 

Since ({r, R)^ is an eigenstate of the momentum/energy operators of the relative and 
CM motions of PI and P2 {not of individual P1/P2), it, evidently, represents a state 
of their mutual superposition (MS). Of course, one can have an alternative picture by 
presuming that each of PI and P2 after their collision at r = falls back on the pre- 
collision side and assumes a kind of self superposition (SS) state {i.e,, the superposition 
of pre- and post-coUision states of one and the same particle). Since PI and P2 on their 
collision exchange their momenta, SS state is also described by ({r, i?)^ and this agrees 
with the fact that PI and P2 are identical particles and we have no means to ascertain 
whether they exchanged their positions or bounced back after exchanging their momenta. 
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Evidently, ({r,R)^ can be used to identically describe the MS state of PI and P2 or 
the SS states of PI or P2 and this observation helps in developing the macro-orbital 
representation for a HC particle in a fluid. 

Analyzing the 1-D motion of two HC particles, we recently, concluded [59] that the 
expectation value of the relative distance between two particles satisfies < x >> A/2. 
Applying this inference to Ck{r)^ state of PI and P2, we find that their < r > should 
satisfy < r >> A/2 for k||r case and < r >> A/2 cos ^ for general case because the 
relative dynamics of two particles, interacting through a central force, is equivalent to 
such motion in 1-D. Thus from the experimental view point two HC particles can reach a 
shortest distance, < r >o= A/2 and in this situation their individual locations (c/. Eqn. 
A-12) are given by < tcmCI) >o— — < i"cm(2) >o— A/4. Using similar result for their 
shortest possible distance in 0— space and < Vncii^) >, etc. we note that Cifc(^)^ state is 
characterized by 

<r> > A/2 or q > g„(= yr/rf), (^-13) 
with d being the nearest neighbor distance of two particles, 

< > > 27r, {A- 14) 

< Vncir) > = < Ad{r) > = 0, {A- 15) 

and 



E{2) = < H{2) > 



+ 



4m 4m 



{A - 16) 



Although, as evident from Eqn. A-16, two HC particles in (^{r,R)^ seems have only 
kinetic energy but this does not mean that HC interaction has no impact on £'(2). We 
note that HC interaction controls £'(2) through q > Qo (Eqn. A-13) and this has been 
demonstrated in Section 2.5 and 2.6. Analyzing the other consequences of Eqn. A-13, we 
note the following : 

(i) Quantum size : A HC particle of momentum q exclusively occupies A/2 space if A/2 > a 
because only then the two particles maintain < r >> A/2. We call A/2 as quantum size 
of a particle. It is evident that quantum size is not an absolute size of a particle, rather 
it depends on the relative momentum k — 2q oi two particles; as such it is the size of one 
particle of momentum q as seen by other particle of momentum —q or vice versa. 

(ii) Zero-point force : It is natural that two particles should experience a repulsive force 
if they try violate < r >> A/2 and this force would come into operation only when our 
system has fixed < r >= d and A/2 is made to increase by some suitable process. In a 
many body system this happens when the system is cooled through a temperature below 
which A/2 tend to cross d or its q tends to fall below its zero-point value qo — n/d. We 
call this force a zero-point repulsion and it can be derived from zero-point energy of the 
particle £<, = h^/8md^ corresponding q — qo- Since this force operates only when the two 
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particles tend have a distance shorter than A/2 and it does not depend of cr, it appears 
that this force represents nothing but the HC repulsion whose operational range gets 
extended from cr to A/2 due to wave nature of particles. 

(iii) Macro-orbital representation : Since two HC particles in ({r, R)^ state either 
experience mutual repulsion (if they somehow have < r >< A/2), or no force (if < r >> 
A/2), it is clear that each of them retains its independent identity. Evidently, each particle 
can be identified to be in its SS state which is represented by a kind of pair waveform 
^ = C='=(r, i?) proposed to be known as macro- or^zta/ described by 

= v^sin[(qi.r,)] exp(Ki.Ri), (A - 17) 

with i {i — 1 or 2) referring to one of the two particles and being equivalent to rcuii) 
(Eqn. A-12) varying from n — Q to ri — A/2, while Ri being the CM position of i— th par- 
ticle. It is evident that each particle in its macro- orbital representation has two motions, 
q and K. The plane wave /T— motion can have any K ranging from to oo, while the 
q is constrained to satisfy Eqn. A-13 due to HC interaction. A macro- orbital identifies 
each particle as an entity of size A/2 moving with momentum K. Since C(r, i?)^ is an 
eigenfunction of the energy /momentum operator of the pair not of individual particle, 
each particle shares the pair energy £^(2) (Eqn. A- 16) equally and we have 

This evidently shows that each particle for its motion behaves as an entity of mass 
4m. 
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Figure -1 : t — T/Tx dependence of different representatives of the order parameter of 
A-transition, viz., nK=o{T*) (Eqn. 37: Curve-A*), nK=o{T) = [I- {T /T^f/'^] (Curve-A), 
n*{T) = [1 -N*{T)/N*{Tx\ (Eqn. 58: Curve-B), order parameter of second order phase 
transition r2(T*) (Eqn. 36: Curve-C*), Q{T) (Eqn. 35: Cnrwe-C), experimental ife — iife 
bond strain {dr — dx)/{do — dx) (Curve-El) and experimental ps{t)/p {Curve-E2). 
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